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Abstract 

The High-Energy Transient Explorer (HETE-2), launched in October 2000, is a 
satellite experiment dedicated to the study of 7-ray bursts in a very wide en- 
ergy range from soft X-ray to 7-ray wavelengths. The intermediate X-ray range 
(2-30 keV) is covered by the Wide-field X-ray Monitor WXM, a coded aperture im- 
ager. In this article, an algorithm for reconstructing the positions of 7-ray bursts is 
described, which is capable of correcting systematic aberrations to approximately V 
throughout the field of view. Functionality and performance of this algorithm have 
been validated using data from Monte Carlo simulations as well as from astrometric 
observations of the X-ray source Scorpius X-l. 
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Introduction 



Although the number of detected 7-ray bursts has increased tremendously 
due to the BATSE instrument on board the Compton 7- Ray Observatory [1], 
unique identifications and quick follow-up observations of the celestial object 
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that harbours the burst site are rare. This issue has been successfully ad- 
dressed by the BeppoSAX satellite, that monitors a large fraction of the sky 
(~ 1 sr) in the intermediate X-ray range [2,3] and is able to derive burst lo- 
calisations with an accuracy of ~ 1'. With such improved localisations, the 
Beppo-SAX team has succeeded in detecting X-ray afterglows [4], by which 
valuable insight into the emission mechanism was gained. Finding 7-ray burst 
counterparts and performing spectroscopy at wavelengths other than the 7-ray 
band during the active phase of a burst as well as providing good localisations 
for follow-up observations is the basic scientific motivation of the HETE-2 
mission. 

This publication is structured as follows: After a brief description of HETE-2's 
instrumentation in section 1, the position reconstruction algorithm is outlined 
in section 2. The performance of the algorithm on Monte Carlo generated 
events as well as on data taken in astrometric observations of the X-ray source 
Scorpius X-l is presented in section 3. In section 4 findings on the point spread 
function are shown. A summary in section 5 concludes the paper. 



1 HETE-2 Mission 

The High Energy Transient Explorer (HETE-2) is a dedicated mission for the 
localisation and spectroscopy of 7-ray bursts. Details of the instrumentation of 
HETE-2 and the mission are given in [5]. Its scientific payload consists of three 
experiments: FREGATE, a scintillation crystal experiment, that is expected to 
provide triggers because of its high sensitivity and large sky coverage, SXC, a 
coded-mask imager based on a X-ray CCD chip design, which is able to localise 
7-ray bursts with very high spatial resolution, and the core experiment WXM. 

The Wide Field X-ray Monitor WXM consists of two perpendicularly oriented 
1-dimensional coded mask cameras, in which photons are detected by position 
sensitive proportional counters. The detectors, one pair for each of the two 
orthogonal systems, are filled with xenon gas at a pressure of 1.4 bar with an 
admixture of 3% carbon dioxide as a quenching gas. Each detector contains 
three carbon wire anodes in which the position of an absorbed photon is in- 
ferred by comparing the accumulated charges at both ends of a wire. WXM 
achieves a position resolution of 1 mm (FWHM) at 6 keV. A veto system be- 
low the counting wires reduces the background due to charged particles. WXM 
is sensitive to photons with energies E in the range 2 keV < E < 30 keV with 
an relative energy resolution of AE/E = 8% at E = 8 keV. 
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The mask pattern, identical for the x- and ^/-system, is placed d = 187 mm 
above the detectors and is composed of 103 elements, each 2 mm in width. On 
third of the elements are open and randomly distributed. The parameters have 
been optimised with respect to source localisation accuracy for the expected 
levels of signal and background photon count rates [6]. The combined x- and 
^/-systems are monitoring a field of view of approximately 60° x 60°. Figure 
1 shows the actual mask pattern of a WXM camera. Further details of the 
WXM cameras can be found in [7a-b] . 

mum ii ■!■ mini hii 

Fig. 1. The coded mask of HETE-2's wide-field X-ray monitor WXM 

In order to provide other experiments with precise localisations of 7-ray burst 
sites within seconds after burst onset, an array of twelve burst alert stations 
has been installed below HETE-2's flight path. Once received, the information 
is relayed to the GRB Coordinate Network (GCN), from where it may be 
obtained by interested observers. 



2 Position Reconstruction Algorithm 

For the sake of readability, the algorithm is described for determining the 
angle of incidence 9 X in the x-detector, completely analogous formulae apply 
to the ^/-system. Figure 2 provides the definition of all quantities involved. 
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Fig. 2. The definition of geometrical quantities in the WXM coordinate system. 



2. 1 Correlation 

The task of localising a single point source is performed by computing the 
correlation function [/ ® m] (x) of the mask pattern m(x) and the recorded 
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intensity distribution f(x); the asterisk denoting complex conjugation: 



The correlation function [/ <g> m]{x) will peak at a value xo, indicating the 
distance f(x) is shifted with respect to m(x). This yields the angle of incidence 
9 X by taxi6 x = xo/d, where d = 187 mm is the distance between mask and 
detector. In m(x), open elements of the mask pattern have been assigned a 
value of 1, whereas closed elements correspond to a value of rj (r — 1) = —1/2, 
where r = 1/3 is the open fraction of HETE-2's imaging system. This method 
is known as balanced correlation and removes the correlation background. 

There exist more sophisticated correlation schemes such as unbiased balanced 
correlation (e.g. [8], [9]), that do take into account partial shadowing of the 
mask onto the detector by bursts in the periphery of the field of view and 
increased coding noise caused by steady strong off-axis sources. For HETE-2's 
flight algorithm we rely on basic balanced correlation for three reasons: The 
width of the mask is almost twice as large as the active area of the detector 
and only bursts at angles > 25° off the optical axis are affected by partial 
shadowing. Furthermore, the running average of the mask's transparency is 
close to constant with deviations being very small, so that the modulation of 
the pattern is insensitive to angle of incidence. Thirdly, in HETE-2's nominal 
survey mode, the satellite is rotated in such a way that strong X-ray sources are 
avoided and the assumption of isotropically incidenting background radiation 
is valid. 



2.2 Imaging Aberrations 

In reality the correlation is complicated by imperfections of the detector: The 
photons travel a finite path before interacting and the loci of photons inter- 
acting in the detector are measured with finite position resolution (in the case 
of WXM to 1 mm, subtending an angle of 18'). Whereas the finite position 
resolution simply results in a blurring of the image, the penetration of ener- 
getic photons effectively induces an additional shift in the image, if the event 
is happening at nonzero angles of incidence, and causes the angles of incidence 
to be overestimated. This effect increases in a complex fashion with distance 
from the center of the field of view and is of the order of 12' at 30° from the 
optical axis for a photon distribution following a power law spectrum with a 
spectral index of V = 1.1. 

In contrast to a time consuming Monte Carlo simulation, in which one deter- 
mines the image of the mask pattern under estimated angles of incidence, the 
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algorithm presented in this article aims at modelling all effects involved in the 
delapidation of the mask pattern, enabling a much faster localisation of the 
7-ray burst site. 



2.3 Convolutions 



Imaging aberrations may be described by applying suitable changes to the 
mask pattern m(x) prior to the calculation of the correlation function. This 
can be archieved by convolving the mask pattern m{x) with integration kernels 
p(x) and r(x), that describe the penetration of highly energetic photons and 
the detector resolution, respectively, 

[p m\ (x) = J m(A) • p(x - X)d\. (2) 



The convolution kernel p(x) that is used for describing the penetration of 
photons is given by equation 3: 



where 9 is the Heaviside function and sgn the signum function The convolution 
yields an altered mask pattern m'{x) = [p0m](i), in which the penetration of 
photons is incorporated. The length scale an of the exponential decay, given by 
equation 4, is equal to the spectrally averaged attenuation length a of photons 
inside the detector, corrected by a projection factor: 

a\\(O x ,e y ) = a. Sfe . (4) 

yl + tan 2 6 X + tan 2 6 y 



Estimates for the angles of incidence 9 X and 9 y follow from a source localisa- 
tion with the unmodified mask pattern m(x) in a first step. The finite spatial 
resolution of the position sensitive proportional counter is described by con- 
volution with a Gaussian integration kernel as in equation 5: 

r(x) = . — exp with c = 1/ V8 ■ In 2 mm. (5) 

v^c I 2c J 



The standard deviation corresponds to a value for full width half maximum 
of 1 mm. In contrast to a, c depends only very weakly on the spectral distri- 
bution of the incident photons and is considered to be constant. After convo- 
lution with r(x), the mask pattern has been modified to m"(x) = [r m'](x). 
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m"(x) is the image of the mask pattern under ideal statistics and the source 
reconstruction with m"(x) should not display any deviation from the ideal be- 
haviour. The correlation [/ ® m"] (x) yields a corrected value 6' x for the angle 
of incidence. 



2.4 Correlation in Fourier Space 



Due to their high numerical complexity it is favourable to perform both the 
convolutions and the image localisation by correlation in the Fourier domain. 
The nomenclature is such that lower case letters denote functions in real space 
and the corresponding upper case letters their Fourier transforms: 

M(k) =J r [m(x)] — — I m(x) exp(— ikx)dx. (6) 

Z7T J 



Convolutions and correlations reduce by virtue of equation 7 to mere multi- 
plications in Fourier space. The asterisk denotes complex conjugation: 

\pOm](x) = ^[P(k)M(k)}, 
\f <g> m"\ (x) = T' 1 [F(k)*M"(k)} . 



Diagram 8 summarises all steps. Starting from the Fourier-transform of the 
mask pattern M(k) = T\m{x)\, both convolutions and the correlation are car- 
ried out in Fourier-space by determining the product C(k) = F(k)*R(k)P(k)M(k). 
Inverse transformation yields the correlation function c(x) = JF~ 1 [C(A;)], from 
which the corrected angle of incidence 9' x is derived: 



/ \ 0p(x) /f \ Qr(x) n( . ®/(x) , \ 

m[x) > m[x) > m [x) >■ c[x) 



(8) 



M(k) > M'(k) > M"{k) > C(k). 

■P(k) -R(k) -F(k)* 



In comparision to the WFC camera on board the BeppoSAX satellite, the 
penetration effect is noticably more pronounced in the case of HETE-2. Eval- 
uating the average attenuation length for a photon at the corner of the field 
of view projected onto the detector plane determines the angular deviation in 
the worst case. This value for WFC is ~ 7 times smaller compared to WXM, 
because of the narrower field of view, the higher gas pressure, which results in 
a shorter average attenuation length and the larger distance between the mask 
and detector. Furthermore, being a genuine 2-dimensional camera, the aber- 
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ration vector field of WFC is radial and does not have the complex angular 
dependence described by equation 4. 



3 Results 

3.1 Monte Carlo Simulations 

3.1.1 Simulated Data 

In order to investigate to which extend the algorithm is capable of correcting 
systematic deviations, a set of 2 14 Monte Carlo images containing a large num- 
ber of photons (n 7 = 3000) was generated, so that the statistical scatter in the 
source reconstruction is small. The simulated burst positions were randomly 
distributed, the angles of incidence ranging from —30° to +30°. The number 
of incident photons was corrected with a factor of \Jl + tan 2 9 X + tan 2 9 y for 
the purpose of coping for the decreasing projected area of the detector at in- 
creasing angles of incidence. The spectral distribution of photons followed a 
power law with spectral index of T = 1.1, which is considered to be shallow, 
but not untypical for 7-ray bursts [10]. Because of the flatness of the spec- 
trum, the aberration effect due to penetrating photons is distinctive and easy 
to study. Localisation of photons by the proportional counter was modelled 
by employing detector response data taken prior to the launch. 

3.1.2 Imaging Aberrations 

Figure 3 shows the deviation 9 X — 0(? ource ) of the uncorrected reconstructed 
burst position from the nominal burst position as a function of angle of inci- 
dence 0( source ), The aberration increases with increasing distance 0( source ) from 
the center of the field of view. Each point corresponds to the localisation 
derived from a Monte Carlo image. 

Furthermore, one observes a widening of the scatter of 9 X — Q(? ource ) with 
increasing distance from the optical axis. This is due to the fact that the aber- 
ration, being proportional to the function a\\(9 x ,9 y ) to first order, is concave 
for 9 X > and convex for 9 X < 0. The magnitude of this curvature increases 
with increasing distance from the origin. 

3.1.3 Performance of the Algorithm 

In order to disentangle the contributions from statistical scatter in the recon- 
struction from residual systematics, the dataset has been split into N = 20 
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distance e (source) from the optical axis [deg] 



Fig. 3. Deviation of the reconstructed event position 9 X from the nominal position 
^(source) ^ & function of the distance 0^ ource "> f r0 m the center of the field of view 
following from a correlation of Monte Carlo generated images with the unmodified 
mask pattern. 

intervals in Q( source ). For each interval i the mean value Si of 9' x — #( source ) as 
well as the standard deviation <jj has been derived. Because the intervals are 
chosen such that they cover an angular range of just 3°, the scatter a, can be 
assumed to be dominated by statistics. An increase in photon number n 7 per 
image would have reduced the statistical scatter, according to <7j ~ 1/ ^fn^. 

From the values Si we derive the quantity S = 1/N ■ \jY.i $h which is the mean 
quadratic deviation from zero. Further on, from <7j we determine the average 
statistical scatter a = 1/N ■ Y^i&i- From those two quantities, we define the 
spot size uo = \J a 2 + S 2 . 

Figure 4 illustrates the performance of the applied corrections: The mean 
quadratic deviation Si as well as the statistical error Uj is plotted against the 
angle of incidence Q(? ource ). As the flatness indicates, Si is consistent with zero, 
i.e. the imaging aberrations have been properly corrected. Moreover, <jj has 
become independent from the angle of incidence #(f ource ) ; which indicates that 
the second order influence of Q( source ) on the reconstruction has been properly 
accounted for. Residual deviations are caused by detector nonlinearities, es- 
pecially in the charge division method by which the position of an absorbed 
photon is determined. 

The mean statistical scatter has been determined to a — 1.53' ± 0.29' and the 
mean systematic deviation from zero to S — 1.00' ± 1.72'. The algorithm thus 
fulfills the design requirements. Figure 5 shows the result of combining the 
reconstructions in x- and y-direction. The scatter plot is circularly symmetric 
and the width of the spot corresponds to uj = 1.83' ± 0.97'. 

In the above analysis, the value for the average attenuation length a has been 
optimised is such a way that S assumed a minimal value. The smallest de- 
viations resulted in the choice a = 3.7 mm, which is in agreement with the 
theoretical value a = 4.0 mm. The theoretical value was derived by determin- 
ing the flux weighted average of tabulated values for the photon attenuation 
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distance e lsou,u,! ' from the optical axis [deg] 



Fig. 4. Deviation of the corrected event position 6' x from the nominal source position 
g(source) &g & f unc ^ on Q f an g\ e Q f incidence Q x source ^ f or the Monte Carlo generated 
dataset. 
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Fig. 5. Scatter plot of all localisations (9' x — 9 x s ° urce \o' y — 9 y source ^) for the Monte 
Carlo set of images. 

length in xenon, where the photon spectrum has been altered by the spectral 
photon absorption probability. 



3.2 Astrometric Observation o/Scorpius X-l 



3.2.1 Scorpius X-l Data 

Data from the X-ray source Scorpius X-l has been taken for calibration pur- 
poses. Due to HETE-2's antisolar pointing, any source will drift through the 
field of view at ecliptic rate, i.e. ~ l°/day. The orientation of the satellite 
was such that the WXM coordinate system was rotated ~ 10° relative to the 
celestial (a, 5) frame. Scorpius X-l was viewed at angles of incidence ranging 
from 9 X = —5° to 9 X = +25° and at 6 y = +5° fixed. The spacecraft attitude 
was known to an accuracy surpassing 1'. 
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As pointed out in [11], the spectrum of Scorpius X-l, being a low-mass X-ray 
binary, can be described by a multicolour blackbody spectrum. Due to the 
integration time of At = 2 sec, short term spectral changes like quasi-periodic 
oscillations average out. Long term variabilities either happen on significantly 
longer time scales compared to the data taking, or they affect energies outside 
HETE-2's sensitivity interval. Determining the average attenuation length by 
a photon flux weighted integral with tabulated values for the attenuation 
length yielded a = 2.5 mm. Thus, the spectrum of Scorpius X-l is much softer 
compared to the previous case. The background subtracted images oconsisted 
roughly of n 7 = 3000 photons. 



3.2.2 Performance of the Algorithm 

The analysis 3.1.3 is repeated for the Scorpius X-l dataset. For each pointing 
5i and cr, have been derived. Figure 6 illustrates the accuracy of the source 
localisation algorithm. Again, the values <5j with their errors cr, are shown as 
a function of angle of incidence Q(? ource ). 
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Fig. 6. Deviation of the corrected event position 6' x from the nominal source position 
^(source) &g & function of angle of incidence Q^ source ^ [ n the WXM-coordinate frame 
for the observation of Scorpius X-l. 

The average quadratic deviation from zero has been determined to 5 = 1.59' ± 
1.65' which is slightly larger than in the previous case, but reflects inaccuracies 
in the satellite pointing. Furthermore, this value is probably underestimated 
because of the incomplete sampling of the field of view. The statistical error 
is marginally smaller due to larger photon numbers and has been derived to 
be a = 1.31' ±0.51'. 

Figure 7 shows the combined localisation in x- and y-direction, converted into 
the celestial coordinate frame. The radius of the spot is u = 2.06' ± 1.31'. The 
average attenuation length a has been varied to yield the minimum value for 
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Fig. 7. Scatter plot of all localisations (a' — a ( source ) ; S' — <j( source )) in celestial coor- 
dinates. Scorpius X-l is located at the origin. 

5. In the case of Scorpius X-l, we find the optimum value to be a = 2.2 mm, 
which corresponds well to the theoretically expected value a = 2.5 mm. 



4 Point Spread Function 



4-1 Autocorrelation Function 



In order to assess the extend to which angular resolving power is affected 
by the source localisation algorithm, the point spread function of the imaging 
system has been derived. For a coded mask camera like WXM, the point spread 
function is equal to the autocorrelation s(x) of the mask pattern m"(x), i.e. 
s(x) = [m" <E> m"\ (x) . To be exact, s(x) describes the isotropised point spread 
function. Anisotropies provoke asymmetries in the central correlation peak of 
the function [m" ® m](x). Again, the determination of s(x) is carried out in 
Fourier-space. Diagram 9 illustrates the derivation of s(x) = T~ x [|M"(/c)| 2 ]: 



/ \ 0p(x) ,/ x 0r(x) „, x ®m"0) , , 

m[x) > m (x) > m (x) >■ s(x) 



(9) 



M(k) > M'(k) > M"(k) ► S(k). 

■P(k) -R(k) | p 



4-2 Variations in Angular Resolving Power 



The width w x of the autocorrelation function s(x) has been determined by fit- 
ting a Gaussian to the central peak. The average width w of the point spread 
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function has been defined as the geometric mean of its extend in x- and y- 
direction (w 2 = w x ■ w y ), because the quantity of interest is the size of the 
spot a source is imaged onto. The value of s(x) at the origin corresponds to 
the correlation strength v x . Accordingly, the average correlation strength v 
has been defined to be equal to the arithmetic mean of v x and v y because of 
its analogy to brightness. 



The angles of incidence have been restricted to 0( source ) = ^je 2 x + Q 2 y < 20°, 
in order to ensure that the mask pattern is fully imaged onto the detector. 
The source positions were arranged on a square lattice with mesh size of 0.5°. 
The average attenuation length a has been set to 3.7 mm, corresponding to a 
power law spectrum with spectral index T = 1.1. 



The left plot in Figure 8 illustrates, that the width w of the autocorrelation 
function increases and how the correlation strength decreases correspondingly 
at larger distances from the optical axis. At increasingly larger distance from 
the optical axis, the functions p(x) and p(y) disperse the correlation peak. 
Because the area underneath the peak has to be conserved, the correlation 
strength v drops accordingly. Thus, the angular resolution power has decreased 



thus by 20% for a burst at 
axis. 



Q(source) 



20° compared to a burst on the optical 



10 



15 



20 



distance |e ,source) | from the optical axis [deg] 



^ 0.95 
> 
sz 

D) 0.90 
© 

u) 0.85 
o 

're 0.80 



5 10 15 20 

distance |9 (soljrce) | from the optical axis [deg] 



Fig. 8. Average width w in terms of FWHM (left) and correlation strength v in 



arbitrary units (right) of the point spread function as a function of distance 
from the optical axis. 



Q{source) 



The scatter in Figure 8 is explained by the fact, that at fixed distance from the 
optical axis one obtains larger values of w for a point situated on a diagonal 
than for a location on one of the axes, because in the latter case only one 
term is affected by the convolution with the functions p(x) and p(y) (formula 
3). An analogous argument applies to the dependence of v on the angles of 
incidence. This effect is of the order of 8% at #( source ) = 20°. 
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5 Summary 

This article describes an alternative correlation algorithm for 7-ray burst local- 
isation with the Vi/XM- device on board the HETE-2 satellite. The performance 
regarding source localisation accuracy has been validated using Monte Carlo 
simulated data and astrometric observations of the X-ray source Scorpius X-l. 

• The algorithm corrects systematic deviations in the imaging process by 
convolving the mask pattern with suitable kernels. The convolution and 
the source reconstruction are both carried out in Fourier space, which is 
favourable with respect to time demand. Currently, the on board recon- 
struction algorithm relies on a library of Monte Carlo generated images. 

• The Fourier algorithm is able to correct the systematic deviation caused 
by photons penetrating into the depth of the detector with high accuracy 
throughout the field of view (compare Table 1). The residual mean quadratic 
deviation is most likely due to nonlinearities in the detector response. 

• The influence of the convolutions on angular resolving power have been in- 
vestigated by determining the point spread function. One observes a widen- 
ing of the correlation peak by ~ 30% together with a drop in correlation 
strength with increasing distance from the optical axis. 

• The algorithm requires but a single, physically meaningful parameter: The 
average attenuation length a. Optimised values for a correspond to theo- 
retically expected values for a source with given spectral properties within 
12% (compare Table 2). Theoretical values follow from formula 10: 

_ JJ<S>{E)p{E,e)a{E)d 2 6dE 
a ~ Jj<5>(E)p{E,6)d 2 6dE 

where $(_E) is the photon spectrum, p(E,9) the probability of detecting a 
photon and a(E) the attenuation length of a photon of energy E. Because 
of the spectral diversity of 7-ray bursts, it is necessary to compute a for 
each event separately. In an on board application the weighted sum over a 
parameterisation of the function a(E) would replace integral 10. 
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(10) 





Monte Carlo data 

oD r VV xllvl 


Scorpius X-1 data 

olJ r VVrllVl 


spot size u 
systematic error 5 
statistical error a 


1.83' ±0.97' 4.31' ±2.28' 
1.00' ±1.72' 2.35' ±4.05' 
1.53' ± 0.29' 3.60' ± 0.68' 


2.06' ±1.31' 4.85' ±3.08' 
1.59' ±1.65' 3.74' ±3.88' 
1.31' ±0.51' 3.08' ±1.20' 



Table 1 



Summary of the localisation performance of the algorithm, stated in standard de- 
viation (SD) and full width half maximum (FWHM), for both the Monte Carlo 
simulation (left column) and the astrometric observation of Scorpius X-1 (right 
column). 





Monte Carlo data 


Scorpius X-1 data 


theoretical 
experimental 


a = 4.0 mm 
a = 3.7 mm 


a = 2.5 mm 
a = 2.2 mm 



Table 2 

Comparison of theoretically expected with optimised values for the spectrally aver- 
aged attenuation length a for different photon distributions, a governs the correction 
for the imaging aberration due to photons penetrating into the depth of the detector. 
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